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The production of mesons in ultra-peripheral collisions of relativistic heavy ions is re-analyzed 
using a projection technique to calculate the amplitudes for the appropriate Feynman diagrams. The 
virtuality of the exchanged photons is fully accounted for in this approach. In the case of two-photon 
fusion, it is explicitly shown that the inclusion of nuclear form factors validates the equivalent photon 
approximation. However, this does not apply to three-photon fusion cross sections. The cross section 
of J/tp production in ultra-peripheral collisions at RHIC and LHC are shown to be much smaller 

. . . than the cross sections for the production of C = even mesons of similar masses. 
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Two-photon physics is the dominant process in e^e^ colliders. This was first shown by Brodsky, Kinoshita and 
Terazawa |l|]. In an earher paper, Low showed that one can relate the particle production by two real photons 
(with energies uji and uj2, respectively) to the particle's decay width, T^j. Since both processes involve the same 
CO ' matrix elements, only the phase-space factors and polarization summations are distinct. Low's formula is 
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C^^ I where M is the particle mass, F-,,^ its decay width, and the delta-function accounts for energy conservation. 

^^ . Another important theoretical development was the realization that the cross sections in colliders are well de- 

L^ ' scribed by replacing the virtual photons by an equivalent field of real photons. One often uses the concept of an 

^^ . equivalent photon number, n{u!), with energy uj. This approximation, called the Weizsacker- Williams method H (or 

r^ ' the equivalent photon approximation) yields for the particle production in colliders ||l| 
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To our knowledge, ref. M was the first to apply a similar approach to study particle production in relativistic heavy 

ion collisions. As compared to e^e~ colliders heavy ions carry the advantage of a larger coupling constant (Za), 

K^ , which increases the cross sections by a large factor. The disadvantage is that one needs to separate the final products 

Vh from those created by strong interaction processes. Inserting eq. |l| in H and using the equivalent photon numbers 

. . 1 appropriate to heavy ions the following expression was obtained in ref. H, to leading logarithmic order, 
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where d = 0.681.., 7 is the Lorentz factor (e.g., 7 = 108 for the RHIC collider at Brookhaven), and i? is a parameter 
which depends on the mass of the produced particle. If M is much smaller than the inverse of a typical nuclear radius, 
then R = 1/M, otherwise R is the nuclear radius. These choices reflect the uncertainty relation in the direction 
transverse to the beam, as explained in ref. H. Since spin 1 particles cannot couple to two real photons H, one 
expects that only spin and spin 2 particles are produced. 

Following these ideas, the two-photon fusion mechanism in heavy ion collisions was exploited by several authors, 
including the possibility to search for the Higgs boson |q-[L||. At present, there are experiments at RHIC/Brookhaven, 
and proposed ones for the Large Hadron Collider at CERN (LHC) llj], which aim to study these phenomena. 
For mesons the cross section is very sensitive to the minimum impact parameter, and refs. 1^,^ have shown that 
corrections to eq. are substantial. These corrections are of geometrical nature and use the equivalent photon 
method, as in eq. |P 

Due to the large theoretical and experimental interest in these phenomena p-H (see also ref. Pq] and references 
therein), it is important to calculate the production mechanism with an alternative approach. We use the projection 



method of ref. [|6| to obtain the meson-production ampHtude in terms of the amphtude for production of quark- 
anti-quark pairs by the time-dependent field of the coUiding nuclei. In section 2 we start with a calculation for 
the production of parapositronium in heavy ions colliders. This will define the calculational steps we need for the 
production of mesons. In particular, we show that the results agree with a recent calculation for this process [[l7[ , 
thus validating the projection method. In section 3 we extend the calculation to the production of C = even mesons. 
In this case, one has to account for the nuclear form factors. We show that the equivalent photon method is obtained 
as a consequence of the cutoff of large photon momenta, imposed by the inclusion of the nuclear form factors. In 
section 4 we calculate the cross section for the production of vector mesons (C = odd) by three-virtual photons. In 
particular, we show that the production rates for J/^ mesons are many orders of magnitude smaller than for the 
C = even mesons of similar masses. 
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FIG. 1. Feynman graphs for two- and three-photon fusion in ultra-peripheral collisions of relativistic heavy ions. 

II. TWO-PHOTON FUSION IN HEAVY ION COLLIDERS 

In the laboratory frame the Fourier components of the classical electromagnetic field at a distance b/2 of nucleus 
1 with charge Ze and velocity (3, is given by (in our notation q = (go, Qt, Qs), and qs = qz) 
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For the field of nucleus 2, moving in the opposite direction, we replace (3 by — /3 and b by — b in the equations above. 
Although /3 ~ 1 in relativistic colliders, it is important to keep them in the key places, as some of their combinations 
will lead to important 7 = (1 — /3^)^^/^ factors. 

The matrix element for the production of positronium is directly obtained from the corresponding matrix element 
for the production of a free pair (see fig. |l|(a)), with the requirement that P+ = P_ = P/2, where P is the momentum 
of the final bound state. That is 



M= M1+M2 

= ~^e^u{ — ) 



d^q 

(2^)4 



^''\^-l) 



'"' m^^.)^A^^H^^.)Pm-.^<-.) 



q2 _ ^,j2/4 



Af2/4 



^(f) 



(5) 



where M is the positronium mass. 

The treatment of bound states in quantum field theory is a very complex subject (for reviews, see [Q). In our case, 
we want to use the matrix element for free-pair production and relate the results for the production of a bound-pair. 
A common trick used in this situation is to convolute the matrix element given above with the bound-state wave 
function. One can show (see, e.g., u%) that this is equivalent to the use of a projection operator of the form 
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where • • • is any matrix operator. The first equation apphes to a spin (parapositronium) and the second to spin 1 
(orthopositronium) particles, respectively. In these equations ^(r) is the bound state wavefunction calculated at the 
origin, and e* is the polarization vector, given by ej_i = (0, l/v^, ±i/\/2, 0) and Cq = (0, 0, 0, 1). 
Using eq. in eq. ra, one gets for the parapositronium production 
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where eA/ii/o- is the antisymmetric Levi-Civita tensor. 

Inserting the explicit form of the electromagnetic fields in eq. R we get 
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The delta functions imply the conditions 

q, = -Po/2/3 and go = -/3P./2 . 
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We also note that Jo^3 = — Js^o a-nd eo^si/ = espOfi ^^nd also that M2, which is obtained by the replacement 
^(1) <-^- A^^^ in A^i, is the same as A^i, i.e., A^i = A^2- In other words, the direct (fig. ||(a)) and the exchange 
Feynman diagrams yield the same result for the matrix element. This is a consequence of the imposed condition that 
P_ = P+ = P/2 and of the projection onto the bound-state. It is an important result that will also show up in the 
diagrams involving three-photons. Gathering all these results, and using £o^3^P,yI^ = |P x I|, we get 
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where E — Pq is the total positronium energy. 

We see that uji and 072 play the role of the (real) photon energies. For real photons one expects 4:UJiLlJ2 = P^, as in 
eq. 0. 

The two-photon fusion cross sections can be obtained by using 
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Since the important impact parameters for the production of the positronium will be 6 > l/rrie ^ R, where R is 
the nuclear radius, the integral over impact parameter can start from 6 = 0. Thus, the integral over impact parameter 
in eq. n5^ yields the delta function 
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We now show that the above equation is equal to the equation obtained in ref . [n7| . First we change the variables 
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It is easy to show that 
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The positronium wavefunction at the origin is very well known. It is given by 
the positronium mass. Thus, eq. |l^ becomes 



da C(3) ctq 



^(0) = M^a^/Mi:, where M is 



(21) 



where 



with 



Jb = / A^(5 (qit + q2t - Pt) dqitrfq2t , 



qit X q2t 

2 2 



A'P 



M^-ql- ql 



C(3) = 1.202. 



and 



CTO 



4ZV 



(22) 



(23) 



We have included the zeta-function ^(3) to take into account the production of the para-positronium in higher orbits, 
besides the production in the K-shell. 

The equation El^ is exactly the same as eq. 2.23 of ref. |1^. Thus we have shown that the approach used in this 
article for the production of a bound particle (in this case, the para-positronium) by means of the two-photon fusion 
yields the same results as in the approach of ref. |jl^ . In that article the total cross section for the production of the 
para-positronium was obtained by separating the regions where a leading order logarithmic approximation could be 
used and a region where the integral in eq. 02 could be solved numerically. To verify their results we will follow a 
different route. Using eq. |l^ we can do the integration over the angle 4> between Pf and qt analytically. We get 
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The triple integral in eq. |2J can be calculated numerically. For RHIC, using 7 = 108 and Au + Au collisions, we find 
a = 19.4 mb. For the LHC , using 7 = 3000 and Pb + Pb collisions, we find a = 116 mb. These are in good agreement 
with the results of ref. 1171 . 



III. PRODUCTION OF C = EVEN MESONS 

We can extend the calculation of the previous section to account for the production of mesons with spin J — Q and 
J = 2 by the two-photon fusion mechanism. The following procedure is to be adopted: 

1. Replace the electron-positron lines by quark-antiquarks in the diagram of figure 1(a). 

2. M in the following formulas will refer to the meson mass. 

3. Replace a^ by a^ (2 J -I- 1) 3^^ Qf, where 3 accounts for the muiiber of colors, and Qi is the fractional quark 
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in terms of T^^, the decay-width of the 



charge. These two last factors will cancel out when we express 

meson. To understand how this is done, lets discuss the basics of the annihilation process of a positronium 
(see also ref. pi|). With probability a^ the e~ can fluctuate and emit a virtual photon with energy We- The 
electron recoils and can travel up to a distance ~ l/rUe (or time ^ me) to meet the positron and annihilate. 
This occurs when e~ and e^ are both found close together in a volume of size (1/me)^, i.e., with a probability 
given by \'^{Q)\^ /m\. Thus, the annihilation probability per unit time (decay width) is F ~ a'^\'^{Q)\^ /m\. 
Angular momentum conservation and CP invariance does not allow the ortho-positronium to decay into an 
even number of photons pi . The description of the annihilation process given above is thus only appropriate 
for the para-positronium. A detailed QED calculation yields an extra 47r in the formula above. This yields 
r^^(^S'o) = 8.03 X 10^ s~^, while the experimental value ||2^ is 7.99(11) x 10^ s^^, in good agreement with the 
theory. For mesons, including the color and the charge factors, as described before, the relationship between 

2 

^(0) and F^^ arise due to the same reasons. One gets F^^ — IGira^ ^(0) /M"^ ■ Sj^iQi- 
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extended to meson decays, should be 



According to these arguments the connection between F^^ and 

valid for large quark masses so that l/rriq ^ \/< r^ >, where v<r^> is the mean size of the meson. Thus, 
it should work well for, e.g. charmonium states, cc. In fact, Appelquist and Politzer [ p4[ have generalized 
this derivation for the hadronic decay of heavy quark states, which besides other phase-space considerations 
amounts in changing a to as, the strong coupling constant. This can be simply viewed as a way to get a 
constraint on the wavefunction |'I'(0)p [g3|. One expects that these arguments are valid to zeroth order in 
quantum chromo dynamics and in addition one should include relativistic corrections. But, as shown in pi| , 
the inclusion of relativistic effects, summing diagrams to higher order in the perturbation series, is equivalent 
to solving the non-relativistic Schrodinger equation. 

4. Change the integration variable to qit and q2t- 

5. Introduce form factors F{qit) and F{q2t) to account for the nuclear dimensions. This is a simple way to eliminate 
the integral over impact parameters and will be justified 'a posteriori', i.e., when we compare our results with 
those from other methods. These form factors will impose a cutoff in qu and q2t, so that 

qit , 92t ^ -^ « M , (27) 

where i? is a typical nuclear size. Taking R = 6.5 fm, we get 1/R ~ 30 MeV. This is much smaller than the 
meson masses. As an outcome of this condition, we can replace Q^ ~ AI^/2 in eq. |l^. 

According the procedures 1-5 we get from eq. n^. 
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We notice that n{uj) is the frequently used form of the equivalent photon number which enters eq. g. Thus, eqs. 
pO] and p^ are the result one expects by using the equivalent photon method, i.e., by using eqs. w and 0. This is 
an important result, since it shows that the projection method to calculate the two-photon production of mesons 
works even for light quark masses (i.e., for tt^). In this case there seems to be no justification for replacing the quark 
masses and momenta by half the meson masses and momenta, as we did for the derivation of eq. g8[ This looks quite 
intriguing, but it is easy to see that the step 3 in our list of procedures adopted is solely dependent on the meson mass, 
not on the quark masses, i.e., if they are constituent, sea quarks, etc. Moreover, the projection method eliminates the 
reference to quark masses in the momentum integrals. The condition p^ finishes the job, by eliminating the photon 
virtualities and yielding the same result one would get with the equivalent photon approximation. 

In the next section we will extend this approach to the calculation of vector meson (J = 1~) production by three 
photons. There we will also apply the results to light quark masses, but we will not be able to check the results 
against the equivalent photon method since we cannot calculate the process as originating from the collisions of three 
real photons. 

We now define a "two- photon equivalent number", J\f2-y{M'^), so that a = (7^^^M2'y{M'^), where 
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To calculate this integral we need the equivalent photon numbers given by eq. ^ The simplest form factor one 
can use for this purpose is the 'sharp-cutoff' model, which assumes that 



F{q^) = 1 for q^ < l/R , and F{q'^) = , otherwise 
In this case, we can use the integral 
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The spectrum possesses a characteristic 1/ui dependence, except for ui 3> 'j/R, when it decreases as l/ui^. 

When the condition 7 '3> MR is met, we can neglect the unity factors inside the logarithm in eq. |34, as well as the 
second term inside brackets. Then, doing the integration of ^ from lj = M'^R/A/^ to w = l/R, we get eq. ^. But, 
eq. pq is an improvement over eq. 0. Eq. is often used in the literature, but it is only valid for 7 ^ MR. This 
relation does not apply to, e.g., the Higgs boson production {Mfjo ^ 100 GeV), as considered in ref. p]. 

For quantitative predictions we should use a more realistic form factor. The Woods-Saxon distribution, with central 
density poj size i?, and diffuseness a gives a good description of the densities of the nuclei involved in the calculation. 
However, this distribution is very well described by the convolution of a hard sphere and an Yukawa function |25|1 . In 
this case, the form factors can be calculated analytically 
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For Au we use R = 6.38 fm, and a = 0.535 fm, with po normalized so that / d^rp{r) = 197. For Pb the appropriate 
numbers are 6.63 fm, 0.549 fm, and 208, respectively [|l9[. With this form factor the two-photon equivalent photon 
number dj\f2-ylduj is also obtainable in a closed form. In table 1 we show the cross sections for the production of 
C = even mesons at RHIC {Au + Au) and LHC {Pb -\- Pb) using the formalism described above. 



Table 1. Cross sections for two- photon production of (C = even) mesons at RHIC {Au + Au) and at LHC (Pb + Pb). 



meson 


mass [MeV] 


r^^ [keV] 


a(+) [nb] 


j^RHIC/^q3 


Mif^/W 


^RHIC [^b] 


a^^^ [mb] 


7I"0 


134 


7.8 X 10-^ 


99 


49 


2.8 


4940 


28 


V 


547 


0.46 


86 


12 


1.8 


1000 


16 


v' 


958 


4.2 


147 


5.1 


1.4 


746 


21 


/2(1270) 


1275 


2.4 


179 


3.0 


1.2 


544 


22 


a2(1320) 


1318 


1.0 


67 


2.9 


1.1 


195 


8.2 


Vc 


2981 


7.5 


8.7 


0.38 


0.7 


3.3 


0.61 


XOc 


3415 


3.3 


2.6 


0.24 


0.63 


0.63 


0.16 


X2c 


3556 


0.8 


2.8 


0.21 


0.56 


0.59 


0.15 



As pointed out in refs. [^jlO|, one can improve the (classical) calculation of the two- photon luminosities by intro- 
ducing a geometrical factor (the 0-function in ref. l9| ) , which affects the angular part of the integration over impact 
parameters. This factor takes care of the position where the meson is produced in the space surrounding the nuclei. In 
our approach the form factors also introduce a geometrical cutoff implying that the mesons cannot be produced inside 
the nuclei. However, it is not easy to compare both approaches directly, as we obtain a momentum representation of 
the amplitudes when we perform the integration over impact parameters to obtain eq. |2^. But we can compare the 
effects of geometry in both cases by using equation Ba. After performing the integral over lu, we can rewrite it as 



cr = / dsC{s)ajj{s) 
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where s — 4:UJiuj2 is the square of the center-of-mass energy of the two photons, cr(s) is given by eq. |lj, and C{s) is the 
"photon-photon luminosity" , given by 
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FIG. 2. Two-photon luminosities (see definition in eq. B7p at RHIC and LHC. Dashed lines include a geometric correction. 



In figure 2 we compare the result obtained by eq. ^ and that of ref . Q . The luminosities for RHIC {Au + Au) and 
for LHC (Pb + Pb) are presented. For RHIC the difference between the two results can reach 10% for very large meson 
masses (e.g. the Higgs), but we notice that for the LHC the two results are practically identical, the difference being 
of the order of 3%, or less, even for the Higgs. Thus, the improved version of eq. ^ given by integrating eq. ^, is 
accurate enough to describe meson production by two- photon fusion. Other effects, like the interference between the 
electromagnetic and the strong interaction production mechanism in grazing collisions, must yield larger corrections 
to the (non-disruptive) meson production cross sections than a more elaborate description of geometrical effects. 

The results in this section are very important for our purpose of calculating the production of vector mesons 
by means of three-photon fusion in peripheral collisions. This could be a relevant process, e.g., for a study of the 
three-photon-vertex in charmonium production. 

One might think that the calculation could be performed by using the equivalent photon approximation that, as 
we have seen in this section, works so well for C = even mesons. However, the introduction of a third photon leads 
to an additional integration, which implies that at least two of the exchanged photons cannot be treated as real ones. 
Nonetheless, the results of this section paves the way to the calculation of production of C = odd mesons. Although 
the use of the projection technique to systems composed of light quarks is questionable, we have seen that it works, 
basically because of the relation G% due to the inclusion of the nuclear form-factors. 



IV. PRODUCTION OF VECTOR MESONS 



Lets now consider the diagram of figure 1(b), appropriate for the fusion of three photons into a C = odd particle. 
According to the Feynman rules, the matrix element for it is given by 
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There will be 12 diagrams like this. But, as we will see below, the upper photon leg in diagram of fig. 1(b) can be 
treated as a real photon, meaning that the equivalent photon approximation is valid for this piece of the diagram. 

One has to use the second of the eqs. ^ to account for the projection onto C = odd particles. The calculation of 
the traces is quite lengthy and was performed using the program FORM |2^]. We have found out that the particle 
is produced with its polarization vector in the transverse direction, as the coefficients accompanying Bq are of higher 
order in I/7. Neglecting such terms we get 
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The above product of the longitudinal components of A^ yields factors proportional to the delta-function, i.e. 
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These delta-functions lead to the conditions 
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where W2 is given by |^. 

The integral over go, ^o and fcs yields a factor 1/2, and the matrix element becomes 

(q* - Pt/2) • e 



(40) 



(41) 



(42) 



Ma = -^{Zay VM*(0) / d-'qtd'-k 



(qt-kt)' + (i?/2/3 + fc3) /7' 



^2 

g2 - M2/4 /j2 _ ^j2/4 
-1 



(P*/2-q,)^+c.2V7^ 



(Pt/2 + kt 



(^3/2 



77^ 



exp (qt • b/2) 



As in eq. O we know that the nuclear form factors imply 



9t , Pt 



R 



<M, 



(43) 



(44) 



(a) 



(b) 



(c) 



i|> 



(d) 



(e) 



# 




(f) 



FIG. 3. Feynman graphs for three-photon fusion in ultra-peripheral collisions of relativistic heavy ions. 

We can thus use the approximation for the propagator {q^ — K'P /4)^^ ~ — 2/M^. Changing the integration variable 
from k to k' = k -I- P/2, the propagators in the second line of ^ become 



(q* + P,/2 - k;)' + {(3EI2 - P3/2 + fc^)' 77^ 



•J2 



'2 , 



V + V/t 



(45) 



The integrand peaks sharply at kj = and one can eliminate k^ from the first term inside brackets. 



Thus, the matrix element WA becomes 



--" - -T '-)' m 



(fqtdkt h 



(qit • e*) exp (q* • b/2) 



X dk:,[kl + klh^ 






(46) 



where qit and (\2t are defined in eqs. |lg and wi and 0^2 are defined in eq. |l4|. 

If this matrix element was the only one being considered it would be easy to show that the one-photon exchange 
with one of the nuclei can be treated in the equivalent photon approximation. This arises from the structure of the 
first term inside the integral in Es. This result was expected in view of our results of the last section. But, the 
two-photon exchange in the lower part of the diagram leads to complicated integrals which cannot be simplified in 
terms of equivalent photons. 

We thus need to calculate the six diagrams which are obtained by interchange of the two-photon lines, as shown in 
figure 3, and multiply the result by 2 to account for the same set of diagrams by inverting the roles of each nucleus. 
We calculate the amplitudes related to the diagrams by means of the same procedures we adopted in eqs. p9^ through 
E^. We find that he matrix elements for the diagrams (b), (c), (e) and (f) of figure 2 are by a factor I/7 smaller than 
those for the diagrams (a) and (d). The amplitude for the diagram (d) in figure 2 yields 



Mo 



? '->' li 



d qtdkt kt 



(W2/W1) (q2t ■ e*) exp (qt • b/2) 



dk-i 



ill + ^!/l^] 
k^ + kl/j^] "^ [qjt + (2cJi -UJ2 + hf/l^] "^ [h +OJ1-UJ2- kl/2uja'^ 



(47) 



The corresponding cross section which is obtained from amplitudes Ma and Aid is given by daa+dab+daint = 2daa- 
The interference term damt yields a contribution of order of 1/7^ after azimuthal integration, and is disregarded. 

The last term in the integrand over k^ dominates the integral in eq. Eq and it is strongly peaked (with width of 
order of Af/7) at k^ ~ uji ~ llJ2- We can replace this value in the other terms of the integrand and take them out of 
the integral. The remaining integral can be done analytically. We get 



Ma = i32 (Za) 



Af3/2 



d^ qtdkt kt 



(qit •e*)exp(qt •b/2) 



1 



1 



kit + ^2/1^] 



[k^ + {lO, L02f/l^] kit + (2^1 - ^lYH 



(48) 



The same trick can be applied to the amplitude of eq. ^. 



We now use eq. |l^ and integrate the squared amplitude over b, multiplying by a factor of 2 to account for the 
amplitude of diagram (d) of figure 3. Again, we insert the nuclear form factors at each of the nuclear vertices to 
account for the nuclear sizes. We also change the integration variables to qit and q2t. The final result, after integrating 

2 

qit ■ e* 



da 
'dK. 



= 10247r 



*(0) {Zaf 



NPE 






dq2t q2t [Fjqjt)]^ 
[ql + (2u;i - c^^r/l^ 



dkt h F{kf) 



{kf + [uj, - oj^Yh^) 



(49) 



We now use the relationship between E and Pz to wi and UJ2 and get rid of the meson wavefunction at the origin. 

2 

The wavefunction \E'(0) cannot be related to the 77 decay widths. But, vector mesons can decay into e^e~ pairs. 

These decay widths are very well known experimentally. Following a similar derivation as for the 77-decay the e+e^ 

decay-width of the vector mesons can be shown ||2^ to be equal to Fg+g- — IGna^ '^'(0) /3A/^ (3 • ^^ Q^). Inserting 

these results in the above equation, the factor (3 • ^^ Q^) will cancel out for the same reason as explained in section 
3, and we get 



da 
duj 



r(-) 



n(w) 



CO 



H{M,u) 



where 



a'^^'^ = 967r 



M3 



with n(cLi) given by pl^ and 



dkt kt F{k?) 



(fc2 + (M2/4CJ - u;)y-/^Y 



(50) 



(51) 



(52) 



The above formulas should also be valid for the production of the ortho-positronium in ultra-peripheral collisions of 
relativistic heavy ions. For RHIC {Au + Au) we obtain a = 11.2 mb, while for the LHC (Pb + Pb) we get cr = 35 mb. 
These numbers are somewhat larger, but are in the same ballpark as the results given in ref. |l^. Note that Coulomb 
corrections [O] substantially modify the positronium production cross section in relativistic heavy ion collisions. This 
is not considered in the present approach. 

In table 2 we present the cross sections for the production of vector mesons by means of the three-photon fusion 
process. We use the form factor given by eq. ^. 

Table 2. Cross sections for three-photon production of vector (C — odd) mesons at RHIC {Au + Au) and at LHC {Pb + Pb). 



meson 


mass [MeV] 


Te+e- [keV] 


a^-^ [nb] 


a^'^J^' [nb] 


a^^^' [nb] 


P' 


770 


6.77 


1740 


137 


1801 


UJ 


782 


0.60 


147 


13 


163 


J/^P 


3097 


5.26 


21 


31 


423 


*' 


3686 


2.12 


5 


12 


155 



We see that the cross sections for the production of vector mesons in ultra-peripheral collisions of relativistic heavy 
ions are small. They do not compare to the production of vector mesons in central collisions. In principle, one would 
expect that the cross sections for three-photon production would scale as {Za)^ , which is an extra Za. factor compared 
to the two-photon fusion cross sections. However, the integral over the additional photon momentum decreases the 
cross section by several orders of magnitude. 
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V. CONCLUSIONS 

We have carried out a derivation of the production of mesons in uhra-peripheral coUisions of relativistic heavy ions 
in terms of a projection procedure. This is useful in order to study the virtuahty content of the exchanged photons. 
We have shown that the cross section for the production of the (para-)positronium is the same as that obtained by 
another calculation p7[ . 

It has also been shown that the inclusion of nuclear form factors leads to cross sections for two-photon fusion which 
agree with those obtained by the equivalent photon approximation. As a byproduct we extended the calculation to 
the production of vector mesons. We have shown that their cross sections are very small and can be neglected for 
practical purposes. 

This work has been authored under Contract No. DE-AC02-98CH10886 with the U.S. Department of En- 
ergy, and partial support from the Brazilian funding agency MCT/FINEP/CNPQ(PRONEX), under contract No. 
41.96.0886.00, is also acknowledged. One of us (CAB) is a fellow of the John Simon Guggenheim Foundation. We 
had also financial support from FAPESP under contract 1999/12987-5. 
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